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ABSTRACT 

Aims. We present a new XVband survey that represents a significant extension to the previous wide-field X^-band imaging survey within the 
0226 - 04 field of the VIMOS-VLT deep survey (VVDS). The new data add ~ 458 arcmin 2 to the previous imaging program, thus allowing us to 
cover a total contiguous area of ~ 600 arcmin 2 within this field. 

Methods. Sources are identified both directly on the final A"-band mosaic image and on the corresponding, deep^ 2 - g'r'f image from the CFHT 
Legacy Survey in order to reduce contamination while ensuring us the compilation of a truly K-selected catalogue down to the completeness limit 
of the X>band. The newly determined A" s -band magnitudes are used in combination with the ancillary multiwavelength data for the determination 
of accurate photometric redshifts. 

Results. The final catalogue totals ~ 52 000 sources, out of which ~ 4400 have a spectroscopic redshift from the VVDS first epoch survey. 
The catalogue is 90% complete down to i^vega = 20.5 mag. We present A^ s -band galaxy counts and angular correlation function measurements 
down to such magnitude limit. Our results are in good agreement with previously published work. We show that the use of A* magnitudes in 
the determination of photometric redshifts significantly lowers the incidence of catastrophic errors. The data presented in this paper are publicly 
available through the CENCOS database. 

Key words, infrared: galaxies - galaxies: general - surveys - cosmology: large-scale structure of Universe 



1. Introduction 

Near-infrared galaxy surveys are widely recognized to give a 
number of advantages with respect to optical surveys as tools 
to study the process of mass assembly at high redshifts. The ob- 
served K-band gives a measure of the rest-frame optical fluxes 
for intermediate redshift galaxies up to z ~ 2, therefore it can 
be more easily related to the galaxy mass in stars. Furthermore, 
the K-band selection leads to the inclusion of extremely red ob- 
jects that would be otherwise missed by a selection in the op- 
tical regime. Additional advantages with respect to optical sur- 
veys are given by the smaller effects of extinction on K-band 

Send offprint requests to: S. Temporin 

* Based on observations collected at the European Southern 
Observatory New Technology Telescope, La Silla, Chile, program 
075.A-0752(A), on data obtained with the European Southern 
Observatory Very Large Telescope, Paranal, Chile, program 070.A- 
9007(A), and on observations obtained with MegaPrime/MegaCam, a 
joint project of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii 
Telescope (CFHT) which is operated by the National Research Council 
(NRC) of Canada, the Institut National des Science de l'Univers of the 
Centre National de la Recherche Scientifique (CNRS) of France, and 
the University of Hawaii. This work is based in part on data products 
produced at TERAPIX and the Canadian Astronomy Data Centre as 
part of the Canada-France-Hawaii Telescope Legacy Survey, a collabo- 
rative project of NRC and CNRS. 



observations and by the smaller required fc-correction, with lit- 
tle dependence on the galaxy types. Near-infrared data are also 
particularly useful for a more accurate determination of photo- 
metric redshifts, a key issue especially in the redshift range 1 
< z < 2, where the measurement of spectroscopic redshifts can 
be challenging. 

In fact, in recent years several efforts have been devoted 
to the compilation of K-selected samples of galaxies, including 
spectroscopic surveys of K-selected sources such as the K20 sur- 
vey ( Cimatti et al.|2002 >. However, surveys reaching very faint 
K-band magnitudes tend to be limited to rather small sky ar- 
eas (and thus are affected by cosmic variance, Labbe et al. (e.g. 
FIRES 2003 down to a depth K ~ 24 over a few square ar- 
cminutes)), while surveys on large areas are often only moder- 
ately deep (e.g. |Daddi et al.|2000||Drory et al. |2001||Kong et al. 
2006, down to K ~ 19 over several hundred square arcminutes). 
A recent overview of the depth and area of published near-IR 
imaging surveys can be found in figure 15 of Forster Schreiber 
[etlLl ( |2rj06l l. 

Only very recently surveys that cover the intermediate 
regime of area and depth have started to appear. Some exam- 
ples are the UKIDSS Ultra Deep Survey and Deep ExtraGalactic 
Survey (DR1; various levels of depth and coverage over sev- 
eral square degrees, | Warren et al. |2007| l, the MUSYC survey 
(5cr^ve ga ~ 21 over 4 fields 10'x 10'each, |Quadri et al.|2007] , 
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and the Palomar Observatory Wide-Fiel d Infrared Survey (1.53 



deg 2 over 4 fields, down to K\/ ega ~ 20.5,|Conselice et al. 



2007). 



Here we present a new /f-band survey that, in the context 



of the VIMOS-VLT deep survey (VVDS |LeFevre et al.|2005), 
enlarges significantly the area already surveyed by |Iovino et al. 
(2005), although to a shallower depth, thus yielding a wide K- 
band contiguous field of ~ 629 arcmin 2 within the 0226 - 04 
VVDS field (F02). This dataset benefits from the multiwave- 
length information already available for the field F02, namely 
optical and ultraviolet imagery available through the VVDS 
(UBVRI; |McCracken eTaLl[2003l |Radovich et aiJ|2004T ) and 
the CFHT Legacy Survey (u*g'r' i'z'), 7 -band imagery available 
for a sub-area of ~ 161 arcmin 2 ( Iovino et al.|2005) , and, in the 
radio regime, 2 .4 GHz VLA and 610 MHz GMRT data ( |Bondi et 
al.|20 03 2007 ). Spectroscopic observations from the first epoch 
VVDS ( |Le Fevre et al.|2005| > targeted 2823 sources from our K- 
band catalogue to A"v ega ^ 20.5. Limiting magnitudes (intended 
as 50% completeness limits) for the available multi-wavelength 
imagery are: u*^ ~ 26.5, g\ K ~ 26.4, r' AB ~ 26.1, t' AR ~ 
25.9, z' ~ 25.0 ( jllbert et al.|2006) ; U AB ~ 25.3 ( |Radovich et 
|al.|2004j, Bar ~ 26.5 V ab ~ 26.2, ^ab ~ 25.9 7ab ~ ^.0 
( |McCracken et al.|2003| ); J AB ~ 24.1 ( |Iovino et al.|2005| ). 

The primary aim of this paper is to describe in detail the 
preparation of our /T-band catalogue and to quantify its reliabil- 
ity and completeness. While the data reduction described in the 
following concerns only the newly obtained, shallower /f-band 
data, the analysis of the properties is carried out on the entire 
/f-band catalogue, which includes the deep part of the survey 
already presented by Iovino et al. ( 2005| ) and can be considered 
complete down to Ky ega < 20.5, as it is shown below. Down to 
this magnitude, the deep part of the survey makes up 26% of the 
catalogue. 

The photometric sample and the spectroscopic sub-sample 
whose general properties are described in this paper are then 
used in a companion paper for the selection and analysis of sam- 
ples of objects with extreme colors ( Temporin et al.|2 008 ). The 



/f-band data described here and in Iovino et al. (2005 ) have sig- 
nificantly contributed to improve the determination of the galaxy 
stellar mass function from the VVDS survey, especially for red- 
shifts z > 1.2 (up to z — 2.5), and for the low-mass tail of 
the function at lower redshifts, z < 0.4 (Pozzetti et "aL~||2007] l. 
Additionally, this /f-band catalogue is well suited to follow the 
evolution of the rest-frame /-band galaxy luminosity function 



(Bolzonel la et al.|[20 08 ). The K-bwd photometry presented in 
this paper has been made available to the astronomical com- 
munity through the CENCOS database ( |Le Brun et al.|2007| ) at 
the URL http://cencosw.oamp.fr/ from where it can be retrieved 
together with the photometry in all other available bands (i.e. 
VVDS photometry in UBVRI(J) and CFHTLS photometry in 
u*g'r' i'z' ■)■ The K-band catalogue we present here can be easily 
obtained by quering the database with the appropriate K magni- 
tude limits. 

2. Observations and data reduction 

New ancillary data to the VVDS 0226-04 field (hereafter F02) 
have been obtained in the K s filter with the SOFI Near Infrared 
Imaging camera (Moo rwood et all] [T998 ) at the ESO New 
Technology Telescope in September 2005 and February 2006. 
These observations cover the region of sky immediately adja- 
cent to that previously targeted by deep J and K s observations 
( Iovino et al.|2005| , with some overlap on the western and south- 
ern sides. Hereafter, we refer to the A^-band simply as Zf-band. 
The observations were done in a series of pointings in a raster 
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Fig. 1. Sky coverage of the /T-band survey. Dotted lines outline 
the pointings nl to n9 of the deep survey ( Iovino et al.||2005] ). 
Dashed lines outline the 24 pointings in the newly surveyed 
area, labelled as in Table [T] The solid line indicates the final 
mosaicked area after excision of the borders. 



configuration, organized in a way to ensure significant overlap 
between adjacent pointings as illust rated in Fig.fT in analogy to 
the observation scheme described in |Iovino'e t al. 12005 ). A total 
of 24 fields ~ 5'x5' in size, with overlapping borders, have been 
observed with a series of jittered 90 s exposures, each obtained 
with a detector integration time DIT=10s, and a number of such 
integrations NDIT=9. Jittering was performed by randomly off- 
setting the telescope within a 30" x 30" box for a total typical 
exposure time of 1 hour per pointing (except for one field which 
was exposed for 1.7 hours) with an average seeingof ~ 171 and 
a pixel scale of 0.288 arcsec pixeL 1 (see Table FT for a list of 
the pointings and the seeing conditions during the observations). 
The airmass of the data ranges from 1.11 to 1.38, except for part 
of one pointing that was observed in February 2006 at an airmass 
~ 1.6. 



After excision of the low signal-to-noise borders and of the 
regions around bright stars (or very nearby extended galaxies) 
the new /^-band images resulted in a newly covered area of 458.2 
arcmin 2 , and, when combined with the previous, deeper obser- 
vations presented in Iovino et al.| ( |2005] ), in a total /T-band area 
of 623 arcmin 2 within F02. Hereafter, we refer to the combined 
deep and shallow A'-band images as to the K-wide image. 

Photometric standard stars from Pers son et al.| ( |1998] ) were 
observed 5 to 9 times per night. Each standard star was cen- 
tered within each quadrant and near the center of the detector 
array through a pre-defined sequence of 5 pointings (DIT=1.2s, 
NDIT=15 each). 

2.1. Data reduction 

The data reduction procedure was largely similar to that detailed 
in |Iovino et al.| ( [2003] >. The individual frames were corrected for 
dark current, flat-fielded, and examined for quality assessment 
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by use of IRAFQ packages and scripts. The quality check in- 
volved the evaluation of fluctuations of the sky-background and 
of the magnitude and the FWHM of a reference star along each 
exposure sequence, as well as a check for the presence of ar- 
tifacts that could compromise the quality of the final coadded 
image. Low-quality images, or images whose seeing was signif- 
icantly worse compared to that of the other images in the same 
sequence were discarded. The final total exposure times and the 
number of useful frames per pointing and observing date, after 
the exclusion of lower quality frames, are listed in Table[T] where 
the individual pointings are named according to the scheme of 
Fig. [T] The average seeing as measured on the final coadded im- 
ages and the photometric correction terms (see below) applied to 
each field are reported, too. 

| |200T| l was 



The IRDR stand-alone software (Sabbey et al 



used for the sky-background subtraction, the coaddition of the 
images that compose each exposure sequence, and the construc- 
tion of the associated weight maps. When the exposure sequence 
of a given pointing was splitted between two different nights, 
two separate coadded images were built, one per night of ob- 
servation. Exposure sequences were splitted during the coaddi- 
tion phase also in the case of a sudden change in the sky con- 
ditions during the observations (e.g. changed background-sky 
level and/or seeing). 

The reduction of the photometric standard stars proceeded 
in a similar way, including the quality assessment of the individ- 
ual frames. A photometric zero-point was determined for each 
night of observation through the measurement of the aperture 
magnitudes of the standard stars observed during the same night 
and their comparison with published magnitudes. Uncertainties 
in the zero-points are of 0.01 - 0.02 mag, depending on the night 
of observation, for the run of September 2005, and 0.03 mag for 
the two observing nights of February 2006. As a value for the 
atmospheric extinction coefficient we adopted the one available 
from the European Southern Observatory web site for the nearest 
dates to our observation^] Thus, the values we used are 0.05 for 
the 2005 observations and 0.03 for the observations in February 
2006. 

The coadded images of F02 were reported to flux values 
above the atmosphere and calibrated with the zero-points that 
compete to their nights of observation. Finally the images were 
rescaled to an arbitrary common zero-point and reported to the 
AB magnitude system with zero point 30.0. Since not all our 
nights were of excellent photometric quality, some refinement 
of the photometric calibration was necessary. To ensure a ho- 
mogeneous calibration across the final mosaicked image, the 
measured magnitudes of the stars in the overlapping regions of 
adjacent pointings and in multiple coadded images of the same 
pointing (obtained either in different nights or with different sky 
conditions), were used to test the photometric calibration and 
to determine the flux scale factors to be applied to the individ- 
ual coadded images, when necessary. For this purpose, we took 
as reference for the photometry the pointing which was observed 
under the best conditions during the observing run, as it emerged 
from our quality checks. This was the pointing SW_21. The zero- 
point uncertainty for that night of observation of SW_21 is of 
0.01 mag. Scaling factors were ignored both in the case of very 
small corrections (< 0.01 mag) and in the case of errors in the 



determination of the scale factor larger than the correction itself. 
No scaling factor was applied to 14 out of 35 coadded frames. 
The correction terms applied to the remaining frames are in the 
range 0.01 - 0.18 mag, with a median value of 0.08 mag; the 
individual values with their uncertainties are given in Table T] 
Finally, after the application of these scale factors, the stars in the 
overlapping regions between the new /T-band observations and 
the previous deep observations were used to check the consis- 
tency of the absolute photometric calibration with that obtained 
for the deep-K survey. No systematic shift was detected in the 
photometric calibration. 

The astrometric calibration of the coadded images was a two- 
step process. A first order astrometric solution was found taking 
as reference the astrometric catalogue of the United States Naval 
Observatory (USNO)-A2.0 ( |Monet |1998] l, then the solution was 
refined by using as reference the position of non-saturated point- 
like sources in the astrometrised CFHTLS /'-band image of the 
field. This procedure, already applied to the de ep-K (jlovino et 
| al.|2005) , BVR ( |McCracken et al.|2003| ), and U ( |Radovich etaL 
2004) images, allowed us to reach a higher relative astrometric 
accuracy and to match sources at the sub-pixel level between op- 
tical and K bands. The rms of these astrometric solutions ranges 
between 07056 and 07123, with a median value of 07099. The 
astrometrised and photometrically calibrated images of the in- 
dividual pointings, weighted by the relevant weight maps, were 
combined in a mosaic together with the deep /f-band images 
and, at the same time, regridded to a final pixel scale of 0.186 
arcsec pixel -1 , in order to match the scale of the CFHTLS opti- 
cal image s. During this op eration, performed with the software 
SWarfP|( Bertin et al.|2002[ ), the relevant flux scale factors, deter- 
mined as explained above, were applied to the individual coad- 
ded images to bring the photometry of the final mosaic to a ho- 
mogeneous basis. A similarly constructed mosaic of the weight 
maps was obtained for use in the source detection phase. See 
|McCracken et al.| ([2003) for further details on the stacking pro- 
cedure with SWarp. The median seeing measured on the K-wide 
image as the FWHM of non-saturated point-like sources is ~ 
170. 



2.1 .1 . Quality of astrometric calibration 

A catalogue of CFHTLS /'-band point-like sources was used to 
verify the goodness of the astrometric calibration of the final K- 
band mosaic. We run SExtractor (Berti n & Arnou ts 1996) on the 
image as a detection algorithm. A representation of the radial 
offsets of K-detected stars with respect to their catalogued posi- 
tions is shown in Fig. [2] No systematic offset has been detected 
either in right ascension or in declination. The measured posi- 
tions agree with the catalogued ones within 0707 (07 12) for 68% 
(90%) of the sources. A separate examination of the residuals 
between Zf-band and /'-band positions as a function of right as- 
cension and declination (Fig. [2j> does not show any dependence 
on the location in the frame. The quality of the astrometric cali- 
bration is uniform throughout the final mosaicked image. The K- 
band positions agree with the /'-band ones within 07047 (07086) 
for 68% (90%) of the point-like sources in both right ascension 
and declination. 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 

2 See SOFI web page: www.ls.eso.org/lasilla/sciops/ntt/sofi/index.html 



3 Software developped by E. Bertin and available at the URL 
http://terapix.iap.fr 
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Table 1. Field F02: Observations summary 



Pointing 


RA 

h m s 


Dec 

o / It 


Obs date 


coadded frames 


texp 

(min.) 


seeing 
(arcsec) 


correction terrr|°] 
(mag) 


S_ll 


02 27 14 


-04 27 03 


6 Sept. 05 


40 


60 


1.11±0.14 


0.032+ 0.014 


S_12 


02 26 57 


-04 27 03 


5 Sept. 05 


40 


60 


0.67+0.09 


0.032+ 0.015 


S_13 


02 26 40 


-04 27 03 


6 Sept. 05 


40 


60 


0.69±0.09 


0.032+ 0.025 


SJ21 


02 27 14 


-04 31 18 


5 Sept. 05 


40 


60 


1.09+0.14 




S_22 


02 26 57 


-04 31 18 


5 Sept. 05 


7 


10.5 


1.11+0.10 




SJ22 


02 26 57 


-04 31 18 


5 Feb. 06 


30 


45 


1.03+0.10 


0.050+0.035 


S_22 


02 26 57 


-04 31 18 


6 Feb. 06 


31 


46.5 


0.96+0.08 


0.050+0.035 


S_23 


02 26 40 


-04 31 18 


6 Sept. 05 


40 


60 


0.79+0.13 


...0 


S_31 


02 27 14 


-04 35 33 


4 Sept. 05 


40 


60 


0.99+0.14 


0.175+0.055 


SJ31 


02 27 14 


-04 35 33 


7 Sept. 05 


10 


15 


0.88+0.06 




S_32 


02 26 57 


-04 35 33 


4 Sept. 05 


39 


58.5 


1.01+0.12 


0.014+0.010 


S_33 


02 26 40 


-04 35 33 


4 Sept. 05 


40 


60 


0.94+0.07 




SW.ll 


02 26 23 


-04 27 03 


7 Sept. 05 


41 


61.5 


0.88+0.09 




SW-12 


02 26 06 


-04 27 03 


9 Sept. 05 


40 


60 


1.09+0.05 




SW_12 


02 26 06 


-04 27 03 


12 Sept. 05 


10 


15 


0.91+0.15 




SW_13 


02 25 49 


-04 27 03 


11 Sept. 05 


20 


30 


1.09+0.01 




SW_13 


02 25 49 


-04 27 03 


12 Sept. 05 


30 


45 


1.14+0.18 




SW_21 


02 26 23 


-04 31 18 


7 Sept. 05 


40 


60 


1.19+0.09 




SW_22 


02 26 06 


-04 31 18 


8 Sept. 05 


40 


60 


1.05+0.15 


0.078+0.021 


SW_23 


02 25 49 


-04 31 18 


12 Sept. 05 


21 


31.5 


1.15+0.16 


-0.142+0.027 


SW_23 


02 25 49 


-04 31 18 


12 Sept. 05 


19 


28.5 


1.06+0.02 


0.078+0.017 


SW_31 


02 26 23 


-04 35 33 


7 Sept. 05 


39 


58.5 


0.81+0.08 




SW_32 


02 26 06 


-04 35 33 


8 Sept. 05 


39 


58.5 


0.85+0.05 


0.102+0.021 


SW_33 


02 25 49 


-04 35 33 


12 Sept. 05 


34 


51 


0.94+0.10 




SW_33 


02 25 49 


-04 35 33 


9 Sept. 05 


10 


15 


1.05+0.05 




W_ll 


02 26 23 


-04 18 33 


6 Sept. 05 


20 


30 


0.85+0.14 


-0.095+0.020 


W_ll 


02 26 23 


-04 18 33 


8 Sept. 05 


20 


30 


1.52+0.15 


-0.116+0.035 


W_12 


02 26 06 


-04 18 33 


9 Sept. 05 


40 


60 


1.40+0.10 


-0.116+0.035 


W_12 


02 26 06 


-04 18 33 


12 Sept. 05 


10 


15 


0.92+0.10 


-0.134+0.040 


W_13 


02 25 49 


-04 18 33 


11 Sept. 05 


40 


60 


1.21+0.10 


-0.116+0.035 


W_13 


02 25 49 


-04 18 33 


12 Sept. 05 


10 


15 


1.10+0.15 


-0.179+0.048 


W_21 


02 26 23 


-04 22 48 


8 Sept. 05 


40 


60 


1.16+0.12 


-0.047+0.020 


W_22 


02 26 06 


-04 22 48 


9 Sept. 05 


40 


60 


1.19+0.07 


-0.078+0.022 


W_22 


02 26 06 


-04 22 48 


12 Sept. 05 


10 


15 


1.18+0.07 


-0.190+0.040 


W_23 


02 25 49 


-04 22 48 


11 Sept. 05 


40 


60 


1.08+0.12 


-0.073+0.049 



" Photometric correction terms that have been applied to bring the photometry into agreement with the reference field SW_21. 



0.4 I i 




36.4 36.5 36.6 36.7 36.8 



Fig. 2. Radial (a), right-ascension (b), and declination (c) residuals between the /T-band and /'-band positions for unsaturated, point- 
like sources. The inner (red) and outer (green) dashed circles/lines enclose 68% and 90% of all objects, respectively. The centroid 
of the distribution of radial residuals is indicated by crossed dashed lines. 
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2.1.2. Quality of photometry 

The quality checks on the image sequences of the target and stan- 
dard star fields described above, and the comparison of the mag- 
nitudes between the deep and the shallow parts of our K-bmd 
survey already gave an indication of the reliability of our pho- 
tometry. 

The photometric errors as a function of the K magnitude 
were estimated by use of repeated exposures on the same region 
of the sky. The source detection was performed with SExtractor 
and the MAG_AUTO parameter from the output file was adopted 
as a measure of the total magnitudes. An example is shown in 
Fig. [3] where two 30 minutes exposures have been used for the 
comparison. Since both images contribute equally and indepen- 
dently to the measured error <x shown in Fig. [3] the uncertainty 
for an individual 30 minutes exposure is given by <r/ y2. Then, 
coadding the two images to obtain the final stacked image (we 
recall that the typical exposure time for the stacked image is 
1 hour) reduces the individual errors by a further factor 1 / y2. 
Hence, the final error will be a half of the value measured with 
the above comparison. Therefore, for the final stacked image of 
a pointing the estimated lex errors are 0.04, 0.12, and 0.17 mag 
at /fvega = 18,19, and 20 mag, respectively. 

As an independent check of the goodness of our photome- 
try, we compared our K-band magnitudes with those from the 
UKIDSS Early Data Release ( |Dye et al. ||2006l> in F02 (Deep 
Extragalactis Survey, DXS, see figure 14 of |Dye et al. (2006) 
for an indication of coverage and depth) and found them to be 
in reasonable agreement, although with some spread. No signif- 
icant systematic offset is present. In Fig. |4]we show the com- 
parison for all cross-identified bona-fide point-like sources that 
were selected for having a SExtractor FLUX JRADIUS param- 
eter r\ji > 2.7 in the /'-band (see Ilbert et aL] [2006 ) and have 
Ky esa < 20 in our catalogue, roughly corresponding to a > 10cr 
detection on our TT-band image. 



3. The A>band catalogue 

3.1. Building the catalogue 

The construction of the catalogue went through a series of steps 
aimed at obtaining a catalogue as complete as possible in the K- 
band - suited to the extraction of /f-selected galaxy samples - 
while minimising the incidence of spurious detections without 
affecting our completeness. This goal was achieved by matching 
the results of two detection procedures, one performed directly 
on the K-bwd image, and the other based on the use of a x 2 
image, constructed using the technique outlined in |Szalay et al. 
( 1999 ). The latter had the main purpose of pruning spurious de- 
tections from the catalogue, especially in the K-fmnl regime, and 
to obtain K-bwd magnitudes already matched with the optical 
ones, which were obtained starting from the same x 2 detection 



image (see also |McCracken et al.||2003| |Radovich et al. 2004 
for the use of x 2 images in the preparation of VVDS photomet- 
ric catalogues in other wavelengths). 

The steps we followed are detailed below. The detection and 
measurement of the sources on the final mosaicked image with 
the associated weight map were obtained with SExtractor by ap- 
plying similar criteria to those described in Iovino et al. (2005). 
As a measure of total magnitudes in our catalogue we use the 
SExtractor parameter MAG_AUTO. 

Step 1. Ax 2 - g'r'i' image built from CFHTLS data has 
been used as the detection image for the measurement of the 
/f-band sources with SExtractor in double-image mode. The use 
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Fig. 3. Photometric error as a function of the /f-band magnitude 
obtained as comparison of the magnitudes from two 30 min ex- 
posures of the same field. The thick solid line indicates the mean 
magnitude difference between the two images and the thin solid 
lines indicate the lcr error as a function of the Ky ega magnitude. 
The typical lcr error in the magnitudes of the final stacked lhr 
exposure will be a half of the error reported in this plot. 
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Fig. 4. Comparison of our photometry with the K-b&nd photom- 
etry from the early data release of the UKIDSS survey ( Dye et al. 



2006). Only point-like sources with a 10cr detection were con- 
sidered. The solid line indicates the mean value of the magnitude 
difference. 



of the x 1 image has the advantage of reducing the number of 
spurious detections, while still picking up most reliable sources 
that are present in our K-bwA image. The x 2 image and the K- 
wide image have sufficiently similar seeing to assure the effec- 
tiveness of this procedure. For an object to be included in our 
catalogue, it must contain at least three contiguous pixels above 
a per-pixel detection threshold of 0.4<x in the^ 2 image. This step 
produced the version 1 (VI) of our catalogue and offered us a di- 
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rect cross-identification of the £"-band sources with those listed 



in the CFHTLS catalogue for the same region of the sky (Ilbert 
|et al.|2006) . The total number 

Step 2. A second version of the catalogue (V2) was produced 
by using SExtractor in single-image mode (i.e. without a sepa- 
rated detection image) on the K-wide image. This was a key step 
for the creation of a Zf-selected catalogue and was necessary to 
ensure that we would not miss any real Zf-band source that es- 
caped detection on the^ 2 image. In this case, to reduce the num- 
ber of fake detections, we chose to include in the catalogue the 
objects containing at least 9 contiguous pixels above a per-pixel 
detection threshold of 1 . 1 cr. 

Step 3. A positional match of sources between the two cat- 
alogues was obtained by adopting l'.'O as an upper limit to the 
separation for a source to be considered the same. Tests with 
different upper limits to the source separation showed that our 
choice was conservative enough to include all sources present in 
both catalogues and sufficiently small to minimize false matches. 
This procedure provided us with a list of objects that were de- 
tected in the K-wide image but were not in the^ 2 image. Besides 
truly red objects, this list included members of (apparent) pairs 
that are resolved in the Zf-band image but appear as an individ- 
ual source in the^ 2 image (see Step 6 below), sources that fall in 
the vicinity of bright objects and are, therefore, masked by bright 
halos in the x 2 image, and, finally, a number of fake detections 
that needed to be expunged. 

Step 4. In order to refine the object list produced at Step 3, 
we evaluated the level of contamination from fake detections as 
a function of magnitude. To this purpose we applied SExtractor 
to the inverse of the K-wide image, with the same parameters we 
used for the production of the catalogue V2. The contamination 
level within 0.5 magnitude bins was estimated as the fraction 
of detections in the inverse image with respect to the original 
image. The resulting contamination rate is 9% in the bin centered 
at Ky esa = 19.75 and reaches 23% in the bin centered at Ky eg - d 
= 20 (i.e. in the range 19.75 - 20.25 mag). We note that this 
contamination rate is relevant for the sources that were detected 
using only the Zf-band image, while it is not representative of 
the contamination from spurious detections of catalogue V 1 , for 
which the effect is minimized by the use of the x 2 image (as 
demonstrated in |Iovino et al.|2005) . 

To limit the amount of contamination, out of all K-hwA de- 
tections for which it was not possible to identify a counterpart 
in the catalogue VI we selected those with a magnitude Kv ega 
< 20.25 and, out of these, we retained only the sources that ap- 
peared reliable upon visual inspection of the K-hwA and optical 
images. These remaining 726 sources were added to the cata- 
logue VI. Their magnitudes in the other available bands were 
determined by running SExtractor on the relevant optical images 
by using the K-hw\d as the detection image. 

Step 5. A positional match of our /f-band detections (cata- 
logue V2) with catalogued VVDS sources in the relevant region 
of the sky, allowed us to ensure the inclusion in our catalogue of 
all objects with a measured spectroscopic redshift. At this step 
we added to our catalogue further 66 K-detected sources. These 
sources did not have a counterpart in the catalogue VI because 
of their vicinity to bright stars in the x 2 image. Also, they were 
missed by our Step 4 because most of them are fainter than the 
threshold ^fvega = 20.25 we adopted for the inclusion of addi- 
tional sources from the catalogue V2. In fact they are ^ 3<x de- 
tections with an average magnitude /^vega ~ 20.9, 13 of them 
having ^fyega < 20.5. 

Step 6. A comparison of the magnitudes of the sources in 
common between the catalogues VI and V2 showed a general 



good agreement, thus indicating that the additional sources iden- 
tified in catalogue V2 could be safely added to the main cat- 
alogue without need of aperture corrections. However, a num- 
ber of sources showed a remarkable magnitude difference be- 
tween the two catalogues and required further investigation. We 
adopted as a criterion to select sources with deviant magni- 
tudes the simultaneous satisfaction of the conditions: |/fv eg a_vi - 
^Vega_v2l > 5<x (where <x was obtained by adding in quadra- 
ture the errors on the individual magnitudes) and |ZTvega_vi - 
^Vega_v2l > 0.2 (to exclude the cases with small absolute devi- 
ations but with very small magnitude errors, that would be se- 
lected by our first condition). The search was limited to objects 
brighter than &/ega = 20.25, in analogy to Step 4. 

About 300 objects were selected this way and visually in- 
spected to check for correct cross-identification between the two 
catalogues and for the presence of possible problems. These ob- 
jects were found to include false matches, matches with /if -band 
false detections, source blends that had been correctly deblended 
only in the A'-band image but not in the x 1 image, and vice- 
versa, source blends that had been correctly deblended in the 
^ 2 -image but remained unresolved as an individual source in the 
A"-band image. We corrected the catalogue for the relevant cases, 
by rejecting false detections and adopting magnitude measure- 
ments based on the A"-band image instead of the ^ 2 -image for 
blends that were unresolved in the latter. Vice-versa, for blends 
that were unresolved in the A"-band but correctly resolved in the 
X 2 -image we kept the identifications and magnitudes of cata- 
logue VI. In both cases the fact that these sources were actually 
blends of multiple sources and not multiple clumps within a sin- 
gle galaxy appeared obvious from the visual inspection of the 
individual optical and near-infrared images. 

Step 7. For the deep area of the survey, a comparison between 
the magn itudes obtained usi ng the new x 2 image and those pre- 
sented in Iovino et al. ( 2005 1, based onaj 2 - BVRIK image, did 
not reveal any systematic difference and showed that the new 
measurements are equivalent (within the errors) to the original 
ones. Therefore, for consistency with the previously published 
paper, in our final catalogue we adopted the original magnitudes 
for the deep part of the survey. 

Step 8. A correction for Galactic extinction on an object-by- 
object basis was applied to the magnitudes in the final catalogue 
by using Schl egel et al.| ( jT998] > dust maps. 



3.2. Completeness 

An indicative estimate of the limiting magnitude of our sur- 
vey can be simply obtained based on the background rms of 
the coadded images and on the seeing during the observa- 
tions. The ncr magnitude limits are given by mag(ncr) = zp - 
2.5 log(ncr VA), where zp is the zero point and A is the area of 
an aperture whose radius is the average FWHM of unsaturated 
point-like sources (see Table [T}. The values estimated this way 
for n = 3, 5 are mag(3cr) ~ 21.4 and mag(5cr) ~ 20.9 (in the 
Vega system). Indeed, by measuring the signal-to-noise ratio as 
a function of K-band magnitude across the final K-wide image in 
a 3H-wide running box with a step of 37"(i.e. 200 pixels) (and 
by imposing the relevant positional constraints to keep the box 
within the borders of the K-wide image as marked with a solid 
line in Fig. [T} we obtained a magnitude limit for 3<x detections 
that varies in the range Ky egR = 21.4 - 22.8, the brightest limit 
being referred to the shallower part of the image, which is the 
subject of the present paper. 
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However, a better characterisation of the photometric prop- 
erties of our final image is given by the estimate of the complete- 
ness level in the source detection at various magnitudes. We de- 
termined the completeness level from our capability to recover 
artificial point-like sources inserted at random positions in our 
image. The detection of artificial sources was performed with 
the same SExtractor parameters adopted for the real sources. 
The procedure we followed does not differ from that adopted 
by Iovino et al.| ( |2005| l and we refer the reader to that paper for 
further details. A representation of the completeness level as a 
function of magnitude for the shallow part of our /T-band survey 
is shown in Fig. [5] In particular we reach a nominal complete- 
ness level of 90% (50%) for point-like sources with Ky esa < 20.5 
(21.5). The completeness level for the deep part of the survey (as 
estimated in |Iovino et al.|2Q05| > is 90% (50%) to K Vega < 20.75 
(22.00). 

However, the determination of the actual completeness of 
our final catalogue is more complex than described above. In 
fact, while the completeness test was run on the /f-band image 
alone, the use of the^- 2 -image played a role in the completeness 
level that was finally achieved. On the other hand, having used 
the x 1 -image for the detection of sources, there is the concern 
that our catalogue is not a purely K-selected one, at least for the 
shallow part of the survey that is the object of this paper. As ex- 
plained in the previous section, for ^vega ^ 20.25 all sources 
with a significant detection in the K-bwd image have been in- 
cluded. Furthermore, for AVega < 20.50 the depth of the images 
composing the x 2 -image implies that the last is fully sufficient to 
recover sources even with the reddest colors while, at the same 
time, reducing the number of spurious detections. Therefore, our 
conclusion is that down to Ky ega - 20.5, where we reach a nom- 
inal completeness level of 90% for point-like sources, we can 
safely state that our catalogue is a truly K-selected one. The 
achieved completeness level is supported by the raw number 
counts (see Sect. 6) that do not fall below a power law at least up 



1 - • 



toK 



Vega 



= 20.5. 



Remaining potential problems for the Ky eg . d > 20.5 regime 
are i) possible cases of association of ^-detected sources with 
correlated noise in the /f-band image and ii) some possible level 
of color incompleteness that could arise from missing extremely 
red, /f-faint sources that could have remained undetected in the 
X 1 image. Therefore, for such faint magnitudes regime we can- 
not state that our catalog is a purely Zf-selected one, having nev- 
ertheless been designed to be as complete as possible in both the 
optical and in K-band. This caveat should be kept in mind when 
in the next sections we refer to our catalogue as a /^-selected 
one: this attribute is, strictly speaking, correct only for Ky e „. d < 
20.5. 



3.3. The K-wide photometric sample 

Our final K-wide catalogue totals 51959 detections. It contains 
22846 objects down to our limiting magnitude (i.e. 50% com- 



pleteness limit) ^Tvega ^ 21.5 and 10615 objects down to our 
90% completeness limit [^J/Tvega ^ 20.5. 

However, to be sure to be dealing with a truly K-selected 
catalogue and to avoid contamination by spurious sources, for 
future analysis we restrict our final photometric sample to the 
8856 objects down to Ky egll < 20.25. The number of catalogued 
objects for different magnitude limits are reported in Table 2]for 
ease of comparison with other existing surveys. At this level, 



4 Hereafter we refer to the 90% and 50% limits as to the "complete- 
ness limit" and the "limiting magnitude" of our survey. 
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Fig. 5. Completeness level of our K-shallow survey as a function 
of magnitude; 50% and 90% completeness levels are marked 
with dashed lines. 

Table 2. Field F02: /^-selected photometric sample 



Aj°^ 19.0 


K< 19.5 


K < 19.75 


K < 20.0 


K < 20.25 


K < 20.5 


3608 


5152 


6180 


6399 


8857 


10615 



" All K magnitudes are expressed in the Vega system. Magnitudes in 
the AB system can be obtained by adding +1.84. 



no distinction has been made between stars and galaxies. This 
issue is addressed later on, in Sect. [5] The magnitude distribution 
of the sample is shown in Fig. [6] The photometric redshifts for 
this sample, determined including the K-band information, are 
presented in Sect. [4] 

3.4. The K-selected spectroscopic sample 

Our shallow Tif-band survey and part of the deep one (i.e. all 
pointings in Fig. [T] except for nl, n2, and n3) are located in the 
four-pass area of the VVDS, wich implies a spectroscopic sam- 
pling rate up to ~ 40% down to 7 A b < 24 ( |Le Fevre et al.|2005| l. 
Our K-wide catalogue (down to the magnitude limit of the sur- 
vey, intended as the 50% completeness limit) contains a total 
of 4059 objects covered by spectroscopic observations, out of 
which 3815 have a successful measure of redshift (81 of them are 
secondary objects with Iab > 24). The actual spectroscopic sam- 
pling rate for the K-wide catalogue down to its 90% complete- 
ness limit is 26.6% (27.9% to K VegR ^ 20.25), while the success 
rate is ~ 95.5% down to the same magnitude limit. Secure (flags 
2, 3, and 4 |Le Fevre et al. 2005) > spectroscopic redshifts are avail- 
able for 1792 galaxies and 23 active galactic nuclei (AGN) down 
to /Tvega ^ 20.25 (12 of these objects have Iab > 24). The num- 
ber of galaxies and stars in the sample within various magnitude 
limits and for various levels of spectral quality is summarised in 
Table [3] The star-galaxy separation for the spectroscopic sample 
relies on the spectral classification. 
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Table 3. Field F02: Number of galaxies/stars in the /^-selected spectroscopic sample 



Flag^] ^ 19.0 K<19.5 K< 19.75 K< 19.8 K<20.0 K<20.25 K<20.5 K<21.5 



any 


893 / 176 


1303/201 


1566/219 


1613/222 


1810/239 


2113/250 


2426 / 265 


3489 / 325 


1,21 


71/7 


126/11 


167/12 


173/12 


199/14 


251 / 14 


312/15 


534/27 


2, 22 


183/11 


304/13 


385/ 14 


394/ 15 


456/20 


551/23 


656 / 27 


1022 / 52 


3,4,23,24 


613/158 


834/ 177 


968/193 


999/ 195 


1097 / 205 


1241/213 


1375/223 


1791/246 


9,29 


7/0 


13/0 


18/0 


18/0 


25/0 


36/0 


48/0 


105/ 


11,211, 12,212 


0/0 


0/0 


0/0 


0/0 


3/0 


3/0 


2/0 


2/0 


13, 14,213,214 


16/0 


21/0 


22/0 


22/0 


22/0 


23/0 


25/0 


26/0 


19,219 


3/0 


5/0 


6/0 


7/0 


7/0 


8/0 


8/0 


9/0 



" Flags indicate the spectral quality and reliability of redshift measurements according to the definitions in Le Fevre et al. (2005 I, 
b All K magnitudes are expressed in the Vega system. Magnitudes in the AB system can be obtained by adding +1.84. In this table the star/galaxy 
separation is based on a purely spectroscopic classification. 
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Fig. 6. Magnitude distribution of our K-wide photometric sam- 
ple down to the magnitude limit of the survey Ky ega < 21 .5 
for all objects (yellow filled histogram) and for objects classi- 
fied as galaxies (red histogram; Sect. 5). The inset shows the K- 
magnitude distribution for the whole Zf-selected spectroscopic 
sample (yellow filled histogram), for the objects classified as 
galaxies within it (red histogram; Sect. 5), and for the galaxy 
subsample with 17.5 < /ab ^ 24.0 (dotted blue histogram). 



The normalized redshift distribution of our /^-selected spec- 
troscopic sample down to the magnitude limits /Tyega ^ 20.25, 
20.50 is shown in Fig. [7] in comparison with the redshift dis- 
tribution of the whole VVDS sample within the same region 
of the sky. The distribution is presented both for the case of 
secure redshifts and for all measured redshifts. We distinguish 
the case in which all available redshifts are considered, irre- 
spective of the /-band magnitudes and the case where only ob- 
jects within the magnitude limits of the spectroscopic survey 
(17.5 < /ab - 24.0) are taken into account. The median red- 
shifts of the distributions down to /Tvega ^ 20.5 closely approach 



0.2 



30.15 



c 

3 0.1 



30.05 


0.2 

30.15 



c 

S 0.1 



; 0.05 




1 1 11 1 1 

any I SB 
flags 1-4, 9 
z_med = 0.78 



z_med 




flags 1-4, 9 
z_med = 0.77 



flags 2, 3, 4 
z_med = 0.76 



z_med 



_L 



17.5 SI AB S24 
flags 2, 3, 4 
z_med = 0.74 



z_med 




Fig. 7. Normalized redshift distribution of the Zf-selected spec- 
troscopic sample (including only sources spectroscopically clas- 
sified as galaxies) down to Ky eg - d < 20.25 (red empty histogram) 
and to /Tvega ^ 20.5 (dashed blue line) compared to the one of 
the entire VVDS sample within the same sky area (cyan filled 
histogram). The left-hand panels include flags 1-4,9, and their 
equivalents, whereas the right-hand panels include only secure 
redshifts with flags 2, 3, 4, 9 and their equivalents. The top pan- 
els show the redshift distribution irrespective of the /-magnitude 
values (i.e. including secondary objects falling outside the se- 
lection criteria of the spectroscopic survey), whereas the bottom 
panels show only the objects within the VVDS magnitude se- 
lection limits. The median redshifts Zmed of the distributions are 
quoted. 



those of the whole VVDS samples both for secure and less se- 
cure redshifts (see values quoted in Fig. [7}. 

The /T-magnitude distribution of the spectroscopic sample is 
shown in Fig. [6] The distributions for both the pure /if-selection 
and the combined K/I ab -selection are displayed. 
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Fig. 8. The / — K vs K magnitude diagram for sources with 
Tab < 24 (black) and 7 A b > 24 (red) and the differential (black) 
and cumulative (green) distributions of K magnitudes illustrate 
the effects of the /ab magnitude cut of the spectroscopic survey 
on the completeness of the Zf-selected spectroscopic sample. A 
slight colour incompleteness arises in the last two 0.5 mag bins 
down to Zfvega ^ 20.5. 



A slight incompleteness in colour is present at faint Zf-band 
magnitudes to the completeness limit of our survey because of 
the /ab ^ 24 magnitude limit of the VVDS specroscopic survey, 
that tends to disfavour faint red objects. The effect of the spec- 
troscopic /ab magnitude limit is illustrated in Fig. [8] where it is 
shown that for /-vega > 19.0 the spectroscopic survey starts miss- 
ing some objects with red I — K colours. Objects with / A b > 24 
represent an increasing fraction of progressively bluer sources 
for increasingly faint K magnitudes. As a result, a ~ 10% (13%) 
colour incompleteness arises in the 0.5 magnitude bin centered 
at vega = 19.75 (20.25). However, photometric redshifts can be 
used to circumvent this problem. 



4. Photometric redshifts 

The photometric redshifts for our Zf-selected sources were ob 
tained with the code Le_Phar^] (developed by S. Arnouts and 
O. Ilbert) , by following the procedure described in Ilbert et 



al. 



2006} and using MAG-AUTO magnitude^] measured in 
the maximum number of available photometric bands in the 
set u'g'r'i'z' BVRIJK. /-band magnitudes (or upper limits) are 
available only for a sub-area of ~ 160 arcmin 2 , i.e. for 26% of the 
objects in our Zf-band catalogue. In all our photometric redshift 
determinations we applied the corrections for systematic offsets 
of the zero-points and the optimised galaxy templates described 
in detail in |Ilbert et aT] ( |2006| l. 



5 http://www.lam.oamp.fr/arnouts/LE_PHARE.html 

6 Magnitudes were measured in all bands with SExtractor in double- 
image mode, by using the^ 2 image as detection image. 



As i t has been sho wn by Ilbe rt et al.| (f2"006), a Bayesian ap- 
proach (Bemtez 2000) with the introduction of an a priori infor- 
mation on the redshift probability distribution function improves 
significantly the determination of photometric redshifts with re- 
spect to the traditional^- 2 method. In this work, consistently with 
Ilbe rt et al.| (j2t)06) and supported by the similarity shown in 
Fig. |7J between K < 20.5 and VVDS redshift distributions, we 
use directly a prior based on the second one and we compare 
the results obtained with/without the application of this prior to 
those obtained with/without the use of the /.-band photometry 
in addition to the other bands during the fitting procedure with 
Le_Phare. This exercise allows us to determine the influence of 
the /.-band on the quality of the photometric redshifts. In com- 
paring the photometric redshifts (zp) with the spectroscopic ones 
(zs), we express the redshift accuracy in terms of cta z /(i +zj ) us- 
ing the normalised median absolute deviation defined as 1.48 x 
median(|Az|/(l + zs)), while for the catastrophic errors we adopt 
the definition \zs - zp\/(l +zs) > 0.15. In analogy to Ilbert et al 



(2006), we splitted our /.-selected spectroscopic sample into an 
/'-bright (17.5 < /' < 22.5) and an /'-faint (22.5 < i' < 24.0) sub- 
sample to check the quality of the photometric redshifts. Also, 
we did our test down to three different limits in K magnitude 
(^Vega i= 20.25, 20.5, 21.5), to evaluate the effects of the inclu- 
sion of /.-faint sources. To ensure the use of reliable photometric 
redshifts, we considered only those based on > 5 photometric 
bands (including the K band). Moreover, to avoid side-effects 
related to the quality of spectroscopic redshifts, we limited the 
comparison to galaxies with spectroscopic flags 3, 4, 23, and 24, 
although we recognize that this choice could bias our subsam- 
ples toward classes of objects whose spectral properties favour 
the redshift measurement. The results are summarized in Table|4] 
where we report the values of cr Az / (I+ , s) and r\ for the selected 
subsamples of galaxies and the total number of objects we used 
for our test. 

Figure [9]illustrates the effects of the inclusion of the /.-band 
and of the prior on the quality of photometric redshifts for the 
Zf-selected spectroscopic sample down to Ky e „. d < 20.25 (the 
results for the samples down to fainter K magnitude limits are 
alike). It is evident that the inclusion of the /.-band, even with- 
out the use of any prior, is more effective than the use of the 
prior alone in reducing the number of catastrophic redshifts, es- 
pecially for objects that are erroneously assigned a photomet- 
ric redshift in the range 2 < zp < 3.5. The best result is ob- 
tained with the use of both the Zf-band and the prior, which gives 
Caz/h+zj) = 0.026 and 77 = 1.61% for the /'-bright subsample, and 
CAz/a+zs) = 0.028 and rj = 2.43% for the entire sample, at least 
in the covered redshift range. 

Objects spectroscopically classified as active galactic nuclei 
are not considered in the comparison shown here. Their inclu- 
sion would increment the catastrophic errors, but would leave 
unchanged the outcome of the comparison, with the best quality 
of photometric redshifts being achieved by including both the 
/.-band and the prior in the fitting process. 

Because of the /ab < 24 magnitude limit of the sample upon 
which the prior is based, we decided to adopt photometric red- 
shifts determined without the prior for sources fainter than this 
limit. In Fig. 10 we show the photometric redshift distribution 
of our sample, after the exclusion of candidate stars, identified 
as described in Sect. [5] The redshift distribution is shown for 
the two magnitude limits /-vega ^ 20.25 and Zfyega ^ 20.5 for 
a purely /.-selected sample and for a sample with an additional 
/-selection that reproduces the magnitude limits of the VVDS 
spectroscopic survey. The effect of the additional /-selection is 
the loss of objects especially at intermediate to high redshift 
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details. The zp — zs line is displayed in green. The results for Le_Phare fits with/without prior and with/without inclusion of /T-band 
photometry are shown with the relevant values of o"a z /(i +ZS ) and 77. The results for the fits with the /T-band are displayed in red. 



17.5 giy, S22. 5 
with K, with prior 

V= 1.69% 




and, therefore, a lowering of the median redshift of the sam- 
ple. The median redshift of the purely Zf-selected photometric 
sample down to K Vsga ^ 20 - 25 ( 20 - 5 ) is zPmed = 0.87 (0.91). 

The photometric redshifts determined in this work are used 
in a companion paper that is dedicated to ^-selected extremely 
red objects (Temp orin et al.||2008 i extracted from our K-bmd 
wide catalogue. 



5. Star-galaxy separation 

We classified point-like objects with a combination of photo- 
metric and morphological criteria, using parameters derived by 
SExtractor and Le_Phare. We adopted a quite complex strategy 
since our aim is to perform the star/galaxy separation even at 
faint magnitudes, where most of the single criteria are known to 
fail. 

The chosen parameters and adopted criteria are the follow- 
ing: 

i) the CLASS_STAR parameter given by SExtractor, provid- 
ing the "stellarity-index" for each object; selected point-like 
sources are obeying CLASS STAR > 0.95,0.90 for objects 
brighter/fainter then i' AB < 22.5 respectively; 

ii) the FLUX_RADIUS parameter, also computed by 
SExstractor in the /T-band and denoted as r\ji, mea- 
suring the radius that encloses 50% of the object total flux; 
we imposed ri/2 < 3.4 pixels as the criterion to be satisfied 
by stellar sources; 

iii) the MILMAX parameter by SExtractor, representing the 
peak surface brightness above the background; the locus of 
selection of stellar objects in the V vs MUJV1AX plot has 
been derived empirically using the spectroscopic sample; 
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Fig. 10. Distribution of photometric redshifts for the ^-selected 
photometric samples down to &/ega < 20.25 (red) and &/ega < 
20.5 (blue). Shaded histograms show the redshift distribution 
when an additional /-selection matching the VVDS selection 
for the spectroscopic survey is applied. The median photomet- 
ric redshifts of each distribution are reported. 



iv) the BzK criterion, proposed by Daddi et al. ( 2004 ), with stars 
characterized by colours z — K < 0.3(5 — z) — 0.5; 
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Table 4. Quality of photometric redshifts determined with/without TT-band photometry and with/without prior down to different 
K magnitude limits and for two intervals of i' magnitudes. For the comparison we used galaxies with good quality spectroscopic 
redshifts (flags 3, 4, 23, 24) and with photometric redshifts derived from at least 5 photometric bands, including the K band. 



K < 20.25 K < 20.5 K < 21.5 K < 20.25 K < 20.5 K < 21.5 
Nobj 890 917 946 347 453 817 



17.5 < (' < 22.5 22.5 < i' < 24.0 



without K, without prior 



o- Az/(1+z) 0.0280 0.0280 0.0280 0.0384 0.0379 0.0351 
tj(%) 7/75 8J8 835 7/78 7/73 7.71 



with K, without prior 



o-Aj/d+a 0.0265 0.0265 0.0267 0.0349 0.0345 0.0327 
??(%) L91 1.96 2.11 5/76 5.30 5.38 



without K, with prior 



o-Az/d+z) 0.0270 0.0268 0.0268 0.0381 0.0367 0.0339 
238 2.62 2.85 61)5 5.52 5.88 



with K, with prior 



o-A Z /(i + z) 0.0264 0.0264 0.0265 0.0350 0.0336 0.0322 
r]{%) L69 1.75 1.90 432 3.97 4.65 



v) the x 2 an d of the SED fitting carried out with Le_Phare 
during the photometric redshift estimate, by using template 
SEDs of both stars (y 2 far )and galaxies (xt a i)> the criterion to 



be satisfied by a stellar source is x 2 



gal 



>0. 



The application of the above criteria to the spectroscopic 
sample is shown in Fig. [TT Objects fulfilling at least 4 of 
the above mentioned criteria have been classified as point-like, 
whereas a less restrictive constraint has been imposed for objects 
brighter than Ky ega = 17. 16, with only 3 over 5 criteria fulfilled. 
Moreover, to be sure not to miss stellar saturated objects, we 
included to the point-like objects selection all the sources with 
r\j2 < 4 and AVega < 14.16 disregarding the other criteria. 

The selection yields 745 point-like objects at Ky eg - d < 20.5, 
with 235 of them belonging to the spectroscopic sample. Within 
this subsample, 80% of the spectroscopic stars are correctly 
identified and only 3 spectroscopic galaxies and 4 active galac- 
tic nuclei with very reliable redshift (flag= 3, 4, 13, 14) fall into 
the star candidate subsample. The galaxy sample is found to be 
contaminated by misclassified stars at the 2% level. 

This method for the star-galaxy separation represents an im- 
provement over both the method used by Ilbert et al. (2006 1 and 
the one adopted by Pozzetti et al. ( 2007 | l, as we verified on the 
spectroscopic samp" 



e. 



Our final classification of the objects takes into account the 
spectroscopic information when this is flagged as highly reliable 
(flag = 3, 4, 23, 24, 13, 14, 213, 214) by overriding the photo- 
metric classification with the spectroscopic one. 



6. Galaxy and star number counts 

Comparing number counts of galaxies and stars with published 
compilations is a good check both of the star-galaxy separation 
efficiency and of the reliability of our photometry, as well as 



the sample reliability and completeness. The differential num- 
ber counts of stars (number 0.5 mag -1 deg~ 2 ) for the F02 wide 
field are shown in Figure 12 To avoid underestimating bright- 
star counts, for this exercise we used the catalogues before ex- 
cising the areas around bright stars. The continuous line is the 
prediction of the model of [Robin et al.| ( [2003| l computed at the 
appropriate galactic latitude. The agreement between observed 
and predicted star counts is very good (the error bars shown are 
Poissonian error bars), confirming the reliability both of our pho- 
tometry and of our star-galaxy separation procedure. 

Figure 13 shows the differential number counts (num- 
ber (0.5 mag) -1 deg~ 2 ) for the F02 wide field. The dotted line is 
obtained by normalising the observed raw number counts to the 
total field after excising areas around bright stars (616 arcmin 2 ). 
The error bars shown are Poissonian error bars and no correc- 
tion for stellar contamination has been applied to the raw counts 
shown. The contamination estimated from the prediction of the 
model of Robi n et al.| ( |2003| l is below 3% for the fainter bins 
shown in the plot. The heavy continuous line shows the total 
galaxy counts obtained after correcting the counts for stellar con- 
tamination estimated using our star-galaxy separation. We did 
not apply any further completeness or contamination correction 
to our data because down to the limiting magnitude (A"v ega = 
20.5) plotted in Figure 13 such corrections are negligible. 



Table [5] lists our raw differential counts, the total raw num- 
ber densities (in units of number (0.5 mag) -1 deg~ 2 ), and the fi- 
nal, corrected for stellar contamination, galaxy densities together 
with their lcr error bars. 

Figure [14] show our total corrected galaxy counts (solid line) 
compared with a selection of literature data. We have followed 
the approach of Cristdbal-Hornillos et al. (2003) and select only 
reliable counts data from the literature, considering only data 
with negligible incompleteness correction and with star-galaxy 
separation applied. We have been conservative in the selection 
of the magnitude intervals plotted in our counts, restricting our- 



12 



S. Temporin et al.: The VVDS: K-band follow-up in the 0226 - 04 field 






o point- like objects 

* spectroscopic stars 

• spectroscopic galaxies 



Fig. 11. The five criteria for the star-galaxy separation (see text) applied to the spectroscopic sample. From top to bottom, and left 
to right, the SExtractor CLA SS-STAR, MU_MAX, and FLUX_RADIUS parameters, the Le_Phare best-fit x 1 values, and the BzK 
diagram (Dad di et al.|2004| are plotted. Each panel illustrates an individual photometric criterion and the corresponding candidate 
stellar sources (red circles), overplotted to spectroscopically classified galaxies (yellow dots) and stars (blue dots). 



selves to bins with relatively large numbers of galaxies, negli- 
gible incompleteness and small contamination corrections. The 
agreement with literature data is very good. Our F02 wide galaxy 
counts confirm the change of slope ar ound Ky^ ~ 18.0 previ- 
ously detected by |Iovino et al.| ( |2005| l, |Gardner et al.| ( |1993| l and 



Cristobal-Hornillos et al. (2003 i. In the range 18 < Ky eg:i < 
20.5 the slope of the galaxy counts is jk ~ 0.30 ± 0.09, while 
in the brighter magnitude range, 15.0 < Ky egll < 18.0, the slope 
is steeper: y K ~ 0.42 ± 0.02. 
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Fig. 12. Differential number counts of stars in the F02 wide field. 
The continuous lines are the prediction of the model of Robin et 

|aT1( |2003l ). 
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Fig. 13. Differential ^-band number counts in the F02 wide 
field. The dotted line shows the raw number densities and er- 
ror bars are Poissonian. The heavy line shows galaxy densities 
obtained after correcting for star contamination. 



7. Clustering analysis for ^-selected data 

In this Section we investigate the clustering properties of point 
like and extended sources in the F02 wide /^-band catalogues. 

We use the projected two-point angular correlation function, 
(l>(6), which measures the excess of pairs separated by an an- 
gle 9,6 + 58 with respect to a random distribution. This statistic 
is useful for our purposes because it is particularly sensitive to 
any residual variations of the magnitude zero-point across our 



Table 5. Column 2 shows the raw differential number counts per 
half magnitude bin. Column 3 shows our raw total number den- 
sities (number 0.5mag~' deg~ 2 ). Column 4 shows our total, dif- 
ferential galaxy densities, corrected for star contamination and 
Column 5 is their poissonian l<x error bar. 
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Fig. 14. The Zf-band galaxy number counts, in units of num 
ber 0.5mag~' deg~ 2 , obtained in this paper are compared to 
a compilation of those from the literature, including counts 
from |Gardner et al.| (|1993|l, |Djorgovski et aTT] <jl995)>, |McLeod| 
et al.| H1995), |Minezaki et al | H1998), |McCracken et al.| d2000)i7 



Kiimmel & Wagner |(2001 1, Saracco et al. (2001 ), Totani et al. 



(|2001),|Maihara et al.| ( |2001|i,|Crist6bal-Hornillos et al.| ( |2003) , 
|Iovino et al.| ( |2005| ), and |Kong et al.| ( |2006| ). 



stacked images. We measure u>(0) using the standard Landy & 



Szalay ( 1993 1 estimator, i.e., 



co e {6) 



DD - 2DR + RR 
RR 



(1) 



with the DD, DR and RR terms referring to the number of 
data-data, data-random and random-random pairs between 9 and 
9 + 69. We use logarithmically spaced bins, with Mog{9) = 0.2, 
and the angles are expressed in degrees, unless stated otherwise. 
DR and RR are obtained by populating the two-dimensional co- 
ordinate space corresponding to the different fields number of 
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Fig. 15. Plot of the correlation function co(9) for the total sample Fig. 16. Results obtained for the amplitude A w at 1 deg of our 
of stars in the /T-band stacks as a function of the logarithm of F02 wide field on galaxy sub-samples of different K-hwA lim- 



the angular pair separation in degrees. At all scales displayed the 
measured correlation values are consistent with zero. 



iting magnitudes. No correction for stellar contamination is ap- 
plied and the error bars on the amplitude are those obtained from 
the fit. 



random points equal to the number of data points, a process 
repeated 1000 times to obtain stable mean values of these two 
quantities. 

7.1. Clustering of point-like sources 

We first measure the angular correlation function co(0) of the stel- 
lar sources. As stars are unclustered, we expect that, if our mag- 
nitude zero-points and detection thresholds are uniform over our 
field, then a>(0) should be zero at all angular scales. 

The results for F02 wide-field ZT-band data are displayed in 
Figure 15 where the correlation function is plotted for the total 
sample of stars obtained from our wide field according to the 
procedure described in section [5] 617 stars in the /T-band inter- 
val 15.00 < K < 20.25. At all scales displayed the measured 
correlation values are consistent with zero. Error bars are ob- 
tained through bootstrap resampling of the star sample (and are 
roughly twice poissonian error bars). 



7.2. Clustering of extended sources 

The procedure followed to measure io{6) is similar to the one 
described above for the star sample. In the case of galaxies a 
positive amplitude of co{6) is expected, and we have to take into 
account the so called "integral constraint" bias. If the real co(6) is 
assumed to be of the form A^O 6 , our estimator ([!]) will be offset 
negatively from the true cl>{9), according to the formula: 



com = Aje- 6 - C) 



(2) 



This bias increases as the area of observation decreases, and 
it is caused by the need to use the observed sample itself to esti- 
mate its mean density, see eg Peebles ( 1980). The negative offset 



AC can be estimated by doubly integrating the assumed true a>(0) 
over the field area Q: 



A/..C 



u)(0)dQ. l dQ. 2 



(3) 



This integral can be solved numerically using randomly dis- 
tributed points for each field: 



C = 



ZN rr (O)0- 
ZN rr (6) 



(4) 



Assuming l"as the pairs minimal scale at which two galaxies 
can be distinguished as separated objects, and 6 = 0.8, we obtain 
the following value: Cwide = 4.50. 

We estimated the amplitude A u for a series of K limited 
galaxy samples by least square fitting A(#~°- 8 -C) to the observed 
cl>(6), weighting each point using bootstrap error bars. Figure 16 
shows the results obtained for 6 different /if -band limiting mag- 
nitudes. No correction for stellar contamination is applied (only 
the objects classified as stars, using the method described in sec- 
tion [5] were excluded from the analysis) and the error bars on 
the amplitude are those from the fit. 

Figure 17 shows the comparison of our results with literature 
data. The continuous, dotted, and dashed lines show the models 
of PLE from Roch e"et"aL"1 ( |1998| l, with scaling from local galaxy 
clustering. 

Table [6] lists our amplitude measurements for the F02 sam- 
ple, in units of 10~ 4 at 1 deg, together with their bootstrap error 
bars. For each limiting magnitude, the total number of objects N 
used in the analysis is also listed. 



8. Conclusions 

We have presented a new Zf-band catalogue that covers a con- 
tiguous sky area of 623 arcmin 2 down to a magnitude limit 
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Fig. 17. We compare our total sample estimate of the amplitude, 
A w at 1 deg, of the angular correlation function with a compila- 
tion of results from the literature in similar magnitude ranges. 



We include mesurements from:|Roche et al. | (| 1 999]», |H>addi et al. 
(p000)l,|Kummel & Wagner 



et al. (2005), and |Kong et aL 



(2001 



Roche et al. } 2003| >7 Iovino 



( 2006 1. The continuous, dotted and 



dashed lines show, for reference, the models of PLE from |Roche| 
et al. ( 1998 ), with scaling from local galaxy clustering. 



Table 6. Observed a>(8) amplitudes A, in units of 10~ 4 at 1 deg, 
together with their lcr error bars, for our F02 wide field. The 
number of objects used in the analysis down to the K magnitude 
faint limits shown in the first column is also listed. 



Magnitude 


0226 


- 04 (F02) 




N 


Al 


dA 


K< 17.25 


561 


40.68 


± 9.24 


K< 17.75 


902 


28.97 


±3.25 


K< 18.25 


1486 


16.92 


±3.32 


K< 18.75 


2427 


17.39 


± 3.41 


K< 19.25 


3661 


13.91 


± 2.67 


K< 19.75 


5470 


13.33 


± 2.66 


K < 20.25 


8059 


13.49 


± 1.57 



K 



Vega 



21.5 and provides us with a 90% complete ^-selected 



sample down to Ky e ga ^ 20.5, although with some possible 
level of color incompleteness potentially affecting extremely red 
sources fainter than Kyega — 20.5, whose inclusion is dis- 
favoured by the procedure to build the catalogue at these faint 
magnitudes (see discussion in Sect. 



of the K band leads to similar advantages as the training on the 
spectroscopic redshift distribution function, namely a consider- 
able reduction of the fraction of catastrophic errors. When no 
a priori spectroscopic information is adopted to train the fitting 
procedure, the additional use of the K band improves signifi- 
cantly the determination of photometric redshifts, as expected. 
The results from this survey show that we expect a significant 
improvement in the accuracy of photometric redshifts obtained 
by future wide-field surveys using near-infrared data (e.g. the 
Visible and Infrared Survey Telescope for Astronomy, VISTArV 

Also, taking advantage of the /T-band photometric parame- 
ters, we have implemented a quite complex set of criteria for 
the star-galaxy separation that improved on methods previously 
adopted by our team (Ilbert et al. 2006 Pozzetti et al. 2007). 
A very good agreement of star and galaxy number counts with 
those present in the literature has proven the effectiveness of our 
star-galaxy separation method as well as the good quality of our 
photometry and the reliability and completeness of our sample. 

The good quality of our photometry as well as the relia- 
bility and completeness of the sample are confirmed by the 
comparison of the number counts of stars and galaxies with 
published compilations. The K-bwA galaxy number counts 
from this work are in excellent agreement with those obtained 
from the ^T-deep sample (Iovino et al. 2005 ) and with a selection 
of compilations from the literature. The projected two-point 
angular correlation function does not show any peculiarity as 
a function of magnitude and angular scale and is broadly in 
agreement with results from the literature. 

The galaxy mass function, the luminosity function, and the 
properties of extremely red objects based on our /T-selected cat 
alogue are presented in companion papers ( Pozzetti et al.|2007 



Bolzonella et al.|2008| |Temporin et al. 2008j |Bondi et al.|2 008). 

A cross-match of this K-bsnd catalogue with radio sources from 



the VVDS-VLA deep field has been carried out bylBondi et al. 



( |2007| l and has revealed a higher incidence of faint /^-band coun- 
terparts (A'vega > 20) among candidate ultra-steep-spectrum ra- 
dio sources with respect to the rest of the radio sources in the 
sample, although based on small number statistics. 
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This is one of the References 



biggest ^-selected samples available to date to this magnitude 
limit and is complemented by UBVRIu*g'r'i'z' photometry - 
available through the VVDS and CFHTLS surveys - as well as 
by VVDS spectroscopy with a sampling rate of ~ 27%, down to 
the 90% completeness limit. 

Good quality photometric redshifts have been obtained for 



the whole sample by following the method outlined by Ilbert et 

[aT" 



2006) and including the /T-band photometry in the fits of 
the galaxy SEDs. By using our spectroscopic subsample we ex- 
plored the effects of the inclusion of /T-band photometry in the 
determination of photometric redshifts. We verified that the use 
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